Pretransitional phenomena in dilute crystals with first-order phase transition 
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Pretransitional phenomena at first-order phase transition in crystals diluted by 'neutral' impu- 
rities (analogue of nonmagnetic atoms in dilute magnets) are considered. It is shown that field 
dependence of order parameter becomes nonanalytical in the stability region of the ordered phase, 
while smeared jumps of thermodynamic parameters and anomalous (non-exponential) relaxation 
appear near transition temperature of pure crystal. 



PACS numbers: 64.60.-i, 05.70.Jk 

In last decades the intensive theoretical and experi- 
mental studies of phase transitions in disordered crystals 
reveal the drastic influence which frozen impurities and 
defects may have on the transition. They can also give 
rise to a variety of pretransitional phenomena. Such phe- 
nomena were first considered by Griffithai in the frame- 
work of dilute Ising model, which has empty sites (with- 
out spins) randomly scattered over lattice. He had shown 
that in such dilute magnet the transition is preceded 
by specific phase, where magnetization is nonanalytical 
function of magnetic field. 

While the bulk transition in this model takes place in 
the (infinite) percolation cluster of spins, the nature of 
pretransitional effects lies in the properties of finite mag- 
netic clusters. Due to the short-range (nearest neighbor) 
interaction between spins all magnetic clusters are com- 
pletely independent and give additive contributions to 
the thermodynamic parameters of the system. Imme- 
diately below the transition point of pure magnet large 
clusters begin to develop high thermodynamic barriers 
between the states with all spins up or down. With due 
reservations these processes can be called 'local phase 
transitions'. As there is a finite probability to find ar- 
bitrarily large magnetic clusters with arbitrarily high 
barriers this results in nonanalytical field dependence 
of magnetizationSi&l. Also same barriers give rise to 
anomalously slow nonexponential relaxation of the order 
parameter in this Griffiths' phased. 

Apparently, the existence of Griffiths' phase is not lim- 
ited to the case of dilute Ising magnet. It should be 
also present in all crystalline solid solutions near generic 
second-order phase transition with short-range order pa- 
rameter interaction, if it takes place in only one compo- 
nent of the solution. The other 'neutral' component (s) 
of the solution (impurities) is supposed to have no influ- 
ence on the transition in the 'pure' copmponent being an 
analogue of nonmagnetic impurities or empty sites in the 
dilute Ising model. Hereafter we refer to this generaliza- 
tion of the notion of dilute magnet as 'dilute crystal with 
a (second-order) transition' or just 'dilute crystal'. 

While the Griffiths' anomalies near second-order phase 
transitions are thoroughly investigated, now nothing is 
known about the pretransitional phenomena near first- 
order transitions in such dilute crystals, in spite of the 
attention paid to a smearing of these transitions by frozen 



disordes2iL§i. But it is evident that some similar anoma- 
lies could exist in this CclSC, clS high barriers between or- 
dered and disordered states of finite clusters of pure com- 
ponent also appear above the transition point of pure 
percolation cluster (if it exists) . 

Let us consider the first-order transition with one- 
component order parameter ip, described in a pure crystal 
by the inequilibrium thermodynamic potential, which is 
even function of <p. We assume that the dimensionless 
(divided by the temperature T) density of such potential 
/ (tp) has a minimum at tp = for T > and two an- 
other minima &ttp ± tp s when T < T + , T_ < T + . Below 
some T , T_ < T < T+, f s = f (±v?s) becomes less than 
fo = / (0), so the first-order transition takes place at To 
in pure crystal. 

It can be shown that in dilute crystal the dependence 
of the order parameter on the conjugate field becomes 
nonanalitycal below T + . To do this we consider the con- 
tributions to the average order parameter from large com- 
pact clusters of pure component with a number of sites 
much greater than that of neighboring impurities and 
diameter much larger than order parameter correlation 
length. Such clusters can be described by the density of 
(inequilibrium) thermodynamic potential of infinite pure 
crystal. Their order parameter distribution in the exter- 
nal field h has the form 

Pn (<p, h) = Z~ l exp [-nf (tp) - nhtp] (1) 



Z n — J dip exp [—nf (tp — nhtp)] 

where n is number of sites in a cluster. It is easy to obtain 
the average order parameter for a cluster 



if)-! 



diptppn (<p,h), 



in a small field 



\h\ < v? s /max(xo,X s ) ■ 



Xo - 1//" (0) , X s = 1//" (<Ps) ■ 
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For n » 1 at T_ < T < T + we have 

X ah + 2<p a (ft) (xs/xo) 1/2 e~ nS f sinh (nh(f s ) 



(2) 
(3) 



l + 2( Xs /xo) 7 e-« 5 /cosh(n^ s ) 

while at T < 1L 

Wi))t,„ = </> s (M- 

Here Sf = f s -f and 

<y9 s (ft) = <£> s tanh (nhtp s ) + \ s h. 



Using the numbers of compact clusters per site, N(n), 
one can find the average order parameter 



n>n 



N(n) (cp), 



Here sum is taken over sufficiently large n > hq for which 
Eqs. (2J and © is valid, 

£max being the largest of correlation lengths of two phases 
(in lattice units) when T_ < T < T + , c > 1 and d is 
space dimension. 

When Sf < the expression in Eq. (JSJ has poles in the 
complex ft plane which are close to that of tanh (nhtp s ) 



± (4m + 1) 



2nip s 



m = 0, 1, .. 



They are arbitrarily close to ft = for large n. As 
N (n) ^ at all n the expansion of average order pa- 
rameter in powers of ft would not converge at T < Tq. 
Above To the expression in Eq. ©has poles at 

2ir 

h rhm ~ ±Sf/(p± i (2m + 1) , m = 0, 1, ... 

nips 

which are close to ±ft c = ±6f/ip s for large n. So at 
T > Tq average order parameter (<p(h)) is nonanalytical 
at the fields ft = ft c and ft — — h c in which transition 
in pure crystal takes place. We can also expect that the 
contributions to (ip (ft)) from sparse clusters would have 
the similar property at larger ft as they have, presum- 
ably, smaller ip s and larger \Sf\ at the same T. Then the 
field dependence of the average order parameter will be 
nonanalytical at all |ft| > ft c . 

The 'local transitions' in finite clusters also manifest 
themselves in the temperature dependence of thermody- 
namic parameters at ft = 0. Consider the compact cluster 
contribution to the zero field susceptibility 



(X)= £ N(n){ X )r,n 



n>n 



(x) T ,n= n {V 2 ) 



(<p)t,, 



At T_ < T < T + we have 

i\T(n) is given by the general formulae 9 , 
tf(n)=5>(n, s)p n (l-p) s 

s 

where 1 — p is the concentration of impurities, and g(n, s) 
is the number of pure clusters with n sites and s neighbor- 
ing impurities. The compactness of clusters considered 
can be accounted for by restriction of the sum to s -c n. 
In the following we will use the lower estimate for N(n), 



N(n) f» p n (1 _p) Smln ("> d ) 



(n,d) = e<jn 



27r d / 2 d d -Vr(d/2) 



(5) 



l/d 



with s m ; n (n, d) being the area of d - dimensional sphere, 
which bounds the volume n. According to Ref. |To|) 
Eq. (J5J describe correctly the main features of large n 
behavior of N(n) either for p — > 1 or p — > except for a 
change of p to some p' > p when p — > (sec Ref. (0)). 
Near Tq, when the condition 



^ e d/(d-l) ^/(d-l) 



n 1 < |<S/| < min jln 

can be fulfilled, we get from Eqs. (jSJ 
^oXo, T > Tq + ST 



1 



1-p 



<x) = 



^OXs 



u 2 ip 2 s , T <Tq-8T 



£ nkN W 

n>n 



Here <5T is defined by the condition no \5f\ ^> 1, so 

<5T ~ T /n 5S, 

where 55 is entropy jump at the transition in pure crys- 
tal. Thus near transition point of pure crystal there is 
slightly smeared jump in the zero field susceptibility. 

Also near To the average entropy of large compact clus- 
ters, 

d 



( S ) = Wr £ N(n)T\nZ n , 



dT 

n>n 

mimics the behavior of pure crystal entropy with prefac- 
tor vq. 

Let us also consider the average modulus of the order 
parameter, (|</?|). This quantity can be interpreted as 
measure of the local order and could be observed in local 
experiments such as NMR. In pure (infinite) crystal it 
becomes nonzero below Tq simultaneously with (ip). In 
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dilute crystal (\tp\) ^ at all T. Indeed, for the compact 
clusters we have 

<M>T,„ « (8xoA«) 1/2 , T > T+ 
(M) T , n ~ fs, T < T_ 

while at T_ < T < T + 

(m)T, n = : — ~ — — —9 — it: ( 6 ) 



l + 2( Xs /xo) 1/2 e- nS f 
Thus in dilute crystal (\<p\) = £ N ( n ) (\<P\) 

n>n 



T,n 



has 



small nonzero value at all temperatures above the tran- 
sition point. 

We cannot describe in the same manner the (additive) 
contributions of the sparse clusters with s > n to the 
thermodynamic parameters. Yet we may suppose that 
they undergo similar 'local transitions', but at tempera- 
tures lower than To due to the lower average coordination 
number. Then continuous sequence of these transitions 
below T will result in steep increase of (\ip\) and decrease 
of (S) at T < To for all concentrations p, while the bulk 
transition takes place at lower temperature (less than 
pTcpi) or vanishes altogether below percolation thresh- 
old. 

We can also show that in the region where phases coex- 
ist large compact clusters have anomalously slow dynam- 
ics with broad spectrum of relaxation times. Here the 
most slow relaxation process is the decay of metastable 
state. In case of one-component nonconserved order pa- 
rameter it normally proceeds via nucleation and growth 
of droplets of stable phas o 12 i 13 . As the nucleation barrier 
must be overcome the time it takes is 

r c = t\ exp A, 



A = max (dl d - x u - l d \6f\) = of- 1 , 

Here a is dimensionless interphase surface tension (di- 
vided by T and multiplied by the area of unit cell face), 
l c is the critical droplet size, 

l c = (d-l)a/\Sf\ 

and T\ is the time inverse proportional to the droplet 
growth ratei£. 

The decay time r c does not depend on system size, 
which is implicitly assumed to be much greater than l c . 
Yet for finite clusters with size less than that of critical 
droplet, n < n c = l d , the above relaxation process trans- 
forms into the sweeping of domain wall through a cluster. 
Then size-dependent barrier, 



A„ = n^/'V - n |<f/|/2, 
should be overcome, which takes the time 
r n = n exp A„. 



(7) 



(8) 



Far from To most clusters exhibit the exponential long- 
time relaxation with single relaxation time r c . When 
T — > T , l c — > 00 and the large part of clusters have size- 
dependent relaxation times. Thus in dilute crystal near 
To the broad spectrum of relaxation times emerges. To 
describe this quantitatively let us consider a contribution 
of compact clusters to a long-time asymptotics of average 
dynamic correlator 



HC(t)) = N(n)C n (t), 

n <n 

C n (t) = (^{t) v (Q)) Tn - (<pf Tn . 

Here C n (t) describes relaxation in n-site compact cluster 
and according to above consideration we have for t — > 00 

c , , _ ( A n exp (-t/r c ) ,n> n c 
™ U \A n exp (-t/r n ) ,n<n c 

Thus near Tq when n c > no 

S(C(t)) = (A) exp (-t/T c ) + B (i) 



B(t)= £ N(n)A 

n 

exp (—t/r n ) 



(9) 



Assuming A n to be a slow varying function of n (slower 
than exp(an)), we can use the steepest descend method 
to find time dependence of B{t) for large t. From 
Eqs. Q-© it follows that at t >> t c clusters with 
= (2/d) d n c , having the largest barrier for a sweeping 
domain wall, give main contribution to B(t), 

B (i) ~ exp (— V T c) j t»T c 



d-l 



When 



Ti exp (a/d) (2l c /d) 



Tq < t « T c , (10) 



-1 ( (d-l)/d\ 

t = ti<t exp [ an J 

the main contribution to B(t) give clusters of size = 
[<T- 1 \n(at/T 1 )] d/ld ~ l) and we get 

B W ~(^) A(P) exp{ln(rt [a-Hniat/n)]^}) 

(11) 

A (p) = -edcr^ 1 In (1-p). 
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Exponential term in falls slower than pure expo- 
nent but faster than all powers of \ jt. Yet for p close to 
1 its variation can be small as compared to that of pref- 
actor in all region IjlUI) and this term should be dropped 
leaving only power-law relaxation. This could happen 
when 



I, 



A(p)<r 
ln(l/j») 



ed- 



La(l-p) 
hip 



>> 1 



Then, not very close to To , for l c << l p we have 
power-law relaxation in the interval l|ll)fl described solely 
by prefactor in J§| , while in the immediate vicinity of T 
for l c >> l p this prefactor can be dropped in (JTTJ. 

Thus near To there appears time interval (|10|) with 
non-exponential relaxation which expands to infinity 
right at To. Owing to the common relaxation mechanism, 
same Eq. Ijllf) holds for long time asymptotic in original 
Griffiths' phase near second-order transition with one- 
component nonconserved order parameter. It was first 
derived in Ref. |5j without power-law prefactor since au- 
thor used the estimate N (n) ~ p n , which does not hold 
for p close to 1. 

Sparse pure clusters with s > n would also exhibit non- 
exponential relaxation of a similar origin. We cannot 
give an exact description of their contributions to the 
dynamic correlator. Yet at every temperature below To 
a 'local transition' takes place in some specific sort of 
sparse pure clusters and their relaxation would dominate 
the infinite time asymptotic. So the above results for 
compact clusters describe the order parameter relaxation 



at t — > oo in close vicinity of To only, while far below it 
relaxation of sparse clusters must be considered. 

Here we may also note that quite similar dynamics 
will appear above the isostructural first-order transitions 
(of the liquid- vapor type), which is described by the po- 
tential with two minima, one at ip — ipo and another 
at <p = (p s , corresponding to high- and low-temperature 
phases. Also the smeared jumps of the susceptibility and 
the entropy appear at To in this case, while the field de- 
pendence of the average order parameter becomes non- 
analytical near To at h c = Sf / (ip s — tpa). 

To conclude, the present qualitative considerations 
show that pretransitional phenomena in dilute crystal 
with first-order phase transition has features rather sim- 
ilar to that of the original Griffiths' phase. They can 
be generalized to the first-order transitions with multi- 
component order parameters and other types of random 
media, which also have cluster size distributions spread- 
ing up to n — oo. So one may expect that analogous 
phenomena could be observed, for example, near melting 
points or the liquid crystal transitions in the substances 
confined in porous media. 
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